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I.  INTRODUCTION 


The  capability  to  predict  and  therefore  understand 
the  discriminatory  features  of  teleseismic  ground  motion 
from  earthquakes  and  nuclear  explosions  is  close  at  hand. 

A  complete  predictive  capability  now  exists  for  estimating 
body  wave  and  surface  wave  teleseismic  ground  motion  caused 
by  an  explosive  source  as  a  function  of  depth  of  burial,  near¬ 
source  environment,  and  earth  structure  (Cherry,  et  al . ,  1973) 
These  advances  in  the  explosion  ground  motion  prediction  area 
are  by  and  large  due  to  a  realistic  calculation  of  the  explo¬ 
sion  equivalent  elastic  source,  based  on  a  specification  of 
the  rock  properties  at  the  source,  and  the  ability  to  propa¬ 
gate  this  source  to  teleseismic  distances  through  a  given 
earth  structure. 

A  similar  predictive  capability  for  earthquakes  will 
exist  when  the  equivalent  elastic  source  for  an  earthquake 
is  determined.  Therefore,  the  central  question  may  be  stated 
as  follows:  Is  it  possible  to  develop  an  earthquake  model, 
whose  parameters  may  be  obtained  from  laboratory  tests  on 
rock  samples,  that  is  capable  of  specifying  the  earthquake’s 
equivalent  elastic  source? 

Our  results  to  date  indicate  that  the  answer  to  this 
question  should  be  affirmative.  The  approach  we  have  taken 
has  involved  incorporating  a  stick-slip  rupture  model  into 
a  two-dimentional  (plane  strain)  Lagrangian  stress  wave 
code.  This  earthquake  model  now  furnishes  the  near  source 
(free-field  and  free  surface)  ground  motion  caused  by  the 
stick-slip  rupture  process. 
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A  description  of  the  model  along  with  free  field 
theoretical  seismograms,  radiation  patterns  and  displacement 
spectra  for  various  fault  lengths,  rupture  velocities  and 
dynamic  stress  drops  have  been  given  by  Cherry  (1973)  and 
we  assume  the  reader  is  familiar  with  that  report.  We  have 
used  the  model  to  simulate  the  February  9,  1971  San  Fernando 
earthquake.  Theoretical  results  for  that  earthquake  are 
presented  in  Section  II  and  compared  with  the  recorded  ground 
motion  at  Pacoima  Dam.  In  Section  III  we  investigate  the 
"far  field"  nature  of  the  displacement  spectrum,  where  in  this 
case  far  field  corresponds  to  ground  motion  five  fault 
dimensions  away  from  the  source. 
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II.  FREE  SURFACE  EARTHQUAKE  GROUND  MOTION 


Two  thrust  fault  calculations  were  run  with  the  only 
variable  being  the  dip  of  the  fault.  The  fault  length  for 
each  calculation  was  5  km  and  the  dynamic  stress  drop  was 
0.5  kbar.  A  homogeneous  geologic  environment  was  assumed 
with  a  P-wave  velocity  of  5.7  km/sec  and  an  S-wave  velocity 
of  3.4  km/sec. 

Figure  1  is  a  schematic  of  the  calculation  (T-45)  in 
which  a  45°  dip  was  assumed  for  the  fault.  The  focal  depth 
was  11.67  km.  The  average  rupture  velocity  was  3.1  km/sec 
as  shown  in  Fig.  2.  The  maximum  relative  displacement  across 
the  fault  was  5.86  meters  and  occurred  4.5  km  from  the  focus 
(Fig.  3). 

Free  surface  ground  motion  was  monitored  at  various 
distances  from  the  epicenter.  Figures  4  through  15  give  the 
vertical  (x)  and  horizontal  (y)  components  of  particle  velocity 

at  the  indicated  epicentral  distances  with  the  origin  of  the 
coordinate  system  at  the  epicenter  (Fig.  1).  In  the  region 
between  the  fault  trace  and  the  epicenter,  the  horizontal 
component  of  particle  velocity  is  always  a  maximum  towards 
the  epicenter  while  the  vertical  component  has  its  largest 
value  in  the  up  direction. 

Figure  16  is  a  schematic  of  the  calculation  (T-90) 
which  assumed  a  90°  dip  for  the  fault.  The  focal  depth  was 
16.5  km.  Both  the  rupture  velocity  and  relative  displacement 
across  the  fault  were  th '  same  as  the  T-45  calculation  and 
are  given  in  Figs.  2  and  3. 

Figures  17  through  28  give  the  vertical  (x)  and 
horizontal  (y)  components  of  particle  velocity  at  the 
indicated  epicentral  distances  for  the  vertical  fault. 
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Free  Surface 

. —  Epicenter 

11.67  km 


Focus 


Fig.  1--Schematic  of  Calculation  T-45.  The 
focal  depth  was  11.67  km  and  t? e 
rupture  length  was  5  km.  The  stations 
monitored  were  at  the  free  surface, 
between  the  epicenter  and  the  fault 
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Fig.  6--Vertical  component  of  particle  velocity  at  x  =  0,  y  =  3.67  km  fo 
Calculation  T-45. 
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Fig.  11 --Vertical  component  of  particle  velocity  at  x  =  0,  y  -  -7.67  kn  for 
Calculation  T-45. 
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Fig.  14- -Horizontal  component  of  particle  velocity  at  x  =  0,  y  =  -11.67  km 
for  Calculation  T-45. 


Calculation  T-45 


Calculation 


jj in , j .. ,  > . .. ' 


20. 


I,  Jf, I,  Li-N  UJA  Iff Jl  11*1 


J..J, ; 


jjni  Jiimi  jifwj.i jwij  ui^wipw^^i 


R-2072 


- ~~~ ♦ 


—  —  ——  ——  ——  —  •  |  NO 


•  • 

•  •  • 

•  O 


•  • 
•  • 
•  • 


L> 

<D 


MMMW  »•••>•  «  "•  "» 


- --4  -- 


_ : - - - -  ° 


(oas/uio) 

XaX3°I9A  9I9XlJ-^d 


23 


. ■*-  - 


»-< 

o 

<u 


e 

,X 


rvi 


O 

II 

X 


nJ 


•H 

L> 

O 


<u 

e 

•H 

H 


0) 

>• 


D 


L) 

•H 

■P 

Jh 

nS 

P. 


m 

o 


4-> 

P  • 

0)  o 
P  CTI 
O  I 
P.H 
£ 

o  P 

u  o 

•  H 
i— I  +-> 

rt  cd 

U  rH 
•H  3 
+->  U 
M  «— I 
s>  rt 
>  u 

i 


o 

CNJ 


GO 

•H 

Ph 


i l 


R-2072 


(oas/uio) 

X^.iooiaA 


24 


Time  (sec) 

Fig.  21- -Horizontal  component  of  particle  velocity  at 
for  Calculation  T-45. 
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Fig.  23- -Horizontal  component  of  particle  velocity  at 
for  Calculation  T-4S. 
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The  maximum  horizontal  component  of  particle  velocity  is 
toward  the  epicenter,  similar  to  the  T-45  calculation. 

However,  the  maximum  vertical  ground  motion  is  down  for  the 
T-90  calculation. 

Figures  29  and  30  give  the  horizontal  and  vertical 
components  of  particle  velocity  recorded  at  Pacoima  Dam 
during  the  February  9,  1971  San  Fernando  earthquake.  The 
maximum  horizontal  component  of  particle  velocity  is  toward 
the  epicenter  while  the  maximum  on  the  vertical  component  is 
in  the  up  direction. 

* 

Estimates  of  the  focal  depth  of  the  San  Fernando 
earthquake  range  between  8-14  km  with  Pacoima  Dam  located 
9-12  km  due  south  of  the  epicenter  (Allen,  et  al . ,  1971, 

Hanks,  1973].  Figures  31  and  32  give  the  horizontal  and 
vertical  components  of  particle  velocity  9.67  km  from  the 
epicenter  for  the  45°  thrust  fault  (calculation  T-45). 

Figures  33  and  34  are  the  ground  motion  components  10  km 
from  the  epicenter  for  the  vertical  thrust  fault  (Calculation 
T-90)  . 

The  character  of  the  ground  motion  at  Pacoima  Dam  is 
best  matched  by  the  T-45  calculation.  Also  the  calculated 
ground  motion  shown  in  Figs.  31  through  34  indicate  that  it 
should  be  possible  to  obtain  a  detailed  match  of  the  first 
few  seconds  of  the  Pacoima  Dam  seismograms  by  increasing  the 
dip  of  the  fault  plane  to  a  value  greater  than  45°. 

It  is  interesting  to  note  that  the  fault  plane  solution 
for  the  San  Fernando  earthquake  gives  a  dip -of  52°  for  the 
fault  (Whitcomb,  et  al . ,  1973)  and  a  maximum  static  displace¬ 
ment  of  5  meters  across  the  fault  near  the  hypocenter. 

(Jangels  and  Frazier,  1973).  Both  results  are  consistent 
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Fig.  30--Vertical  component  of  particle  velocity  at 

Pacoima  Dam  from  the  San  Fernando  earthquake 
Motion  in  the  positive  direction  corresponds 
to  ground  velocity  up. 
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32--Vertical  component  of  particle  velocity  9.67  km  from  the 
epicenter  from  Calculation  T-45.  Positive  direction 
corresponds  to  ground  velocity  down. 
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Fig.  34--Vertical  component  o£  particle  velocity  10  km  from  the 
epicenter  from  Calculation  T-90.  Positive  direction  is 
ground  velocity  down. 
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with  our  calculations.  Also,  we  agree  with  Hanks'  (1973) 
conclusion  that  the  San  Fernando  earthquake  was  initiated 
with  massive  but  localized  rupture  in  the  hypocentral  region. 
There  seems  little  doubt  that  the  rupture  is  capable  of 
being  simulated  by  the  stick-slip  process  (Cherry,  1973) 
with  the  parameters  in  the  model  specified  by  data  from 
laboratory  tests  on  rock  samples. 
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III.  FREE  FIELD  DISPLACEMENT  SPECTRA 


The  free  field  displacement  spectra  reported  previously 
(Cherry,  1973)*  were  obtained  at  distances  of  the  order  of 
two  fault  dimensions  from  the  center  of  the  fault  and  were 
contaminated  by  near  field  effects.  In  an  attempt  to  uncover 
the  nature  of  the  free  field  displacement  spectrum  we  performed 
a  calculation  and  monitored  free  field  ground  motion  five  fault 
dimensions  away  from  the  center  of  the  fault. 

Figure  35  is  a  schematic  of  the  calculation.  The  fault 
length  was  1  km.  The  dynamic  stress  drop  was  0.25  kbar  and 
the  rupture  velocity  equaled  2.3  km/sec.  The  maximum  relative 
displacement  across  the  fault  was  0.247  m  and  occurred  0.35  km 
from  the  focus. 

Figure  36  is  the  spectrum  of  the  y  component  of 
displacement  at  R  =  5  km,  0  =  180°.  The  high  frequency  portion 
of  the  spectrum  decays  as  f  3 .  The  cornel  frequency  is 
approximately  1.9  Hz  (10  ,Z8)  and  the  spectrum  is  flat  in  the 
interval  0.2  Hz  <  f  <  1  Hz.  The  relation  established  by 
Cherry  (1973)  between  corner  frequency  (fc),  fault  length  (L) 


Tt - 

In  that  report  the  FFT  of  particle  velocity  was  divided  by 
frequency  in  Hz  to  obtain  the  displacement  spectrum. 

This  division  introduces  a  factor  of  2tt  in  the  spectrum. 
Therefore  if  u(t)  is  ground  displacement,  then  the  dis¬ 
placement  spectrum  u  is  given  by 

00 

U  =  2-n J"  u(t)e"ia)t  dt 
-  00 

This  definition  of  the  displacement  spectrum  is  maintained  in 
this  report. 
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and  rupture  velocity  (VR)  at  this  azimuth  is 

fc  =  0.16  *  0.7  (JJfi)  . 

This  gives  fc  =  1.77  Hz  for  Vr  »  2.3  km/sec  and  L=  1  km. 

Figure  37  is  the  displacement  spectrum  for  the  x 

component  of  displacement  at  R  =  5  km,  0  =  100°.  The  spectrum 
is  similar  to  that  of  Fig.  36  except  the  corner  frequency 
has  shifted  to  0.9  Hz  CIO'0- 0S)  due  to  a  Doppler  effect  ’ 
from  the  finite  rupture  velocity.  At  this  azimuth  Cherry 
(1973)  obtained  the  relation 

fc  =  0.15  ♦  0.3  (^)  . 

This  gives  fc  =  0.84  Hz  for  VR  =  2.3  km/sec  and  L  -  1  km. 

Figure  38  is  the  spectrum  at  the  epicenter  for  the 

vertical  fault  discussed  in  Section  II  (Calculation  T-90) 
and  shows  the  effect  of  the  free  surface  in  shaping  the 
spectrum.  Apparently  spectral  peaking  in  the  epicentral 
region  is  caused  by  Rayleigh  wave  generation  at  the  free 
surface,  which  provides  a  mechanism  for  reducing  the  static 
displacement  from  twice  its  free  field  value.  Since  the 
high  frequency  components  of  the  incident  wave  are  amplified 
by  the  free  surface  a  peaked  spectrum  results. 

Therefore,  we  should  be  careful  in  deducing  the  nature 
of  the  free  field  spectrum  from  near  field,  surface  measure¬ 
ments  until  the  effect  of  the  free  surface  on  the  spectrum 
is  understood.  This  is  especially  true  for  thrust  and  normal 
faults  where  the  incident  wave  contains  large  P  and  SV 
components.  The  SH  component  of  ground  motion  at  Bear 
Valley  should  be  heavily  weighted  when  the  spectral  content 
of  the  source  is  deduced  from  field  data. 
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Fie.  38- -Displacement  spectrum  for  the  horizontal  component  of  ground 
motion  at  the  epicenter  for  Calculation  T-90.  The  station 
was  2.8  fault  dimensions  from  the  center  of  the  fault. 
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